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ABSTRACT 

Observations of planetary nebulae (PNe) by Sterling, Dinerstein and Bowers have revealed abundances in 
the neutron-capture element Germanium (Ge) from solar to factors of 3 - 10 above solar. The enhanced Ge is 
an indication that the s/ow-neutron capture process (s process) operated in the parent star during the thermally- 
pulsing asymptotic giant branch (TP-AGB) phase. We compute the detailed nucleosynthesis of a series of AGB 
models to estimate the surface enrichment of Ge near the end of the AGB. A partial mixing zone of constant 
mass is included at the deepest extent of each dredge-up episode, resulting in the formation of a 13 C pocket in 
the top ~ 1 /10 th of the He-rich intershell. All of the models show surface increases of [Ge/Fe] < 0.5, except the 
2.5M Q , Z = 0.004 case which produced a factor of 6 enhancement of Ge. Near the tip of the TP-AGB, a couple 
of extra TPs could occur to account for the composition of the most Ge-enriched PNe. Uncertainties in the 
theoretical modeling of AGB stellar evolution might account for larger Ge enhancements than we predict here. 
Alternatively, a possible solution could be provided by the occurrence of a late TP during the post- AGB phase. 
Difficulties related to spectroscopic abundance estimates also need to be taken into consideration. Further study 
is required to better assess how the model uncertainties affect the predictions and, consequently, if a late TP 
should be invoked. 

Subject headings: stars: AGB and post-AGB stars — planetary nebulae: general — nuclear reactions, nucle- 
osynthesis, abundances 



1. INTRODUCTION 

After the TP-AGB phase is terminated, low to intermediate 
mass stars (with initial masses ~ 0.8 to 8M Q ) evolve through 
the PNe phase before ending their lives as white dwarfs (see 
the recent review by Herwig 2005, and references therein). 
The whole envelope is lost by low-velocity stellar winds dur- 
ing the AGB phase. The gaseous nebula is the remnant of 
the deep convective envelope that once surrounded the core, 
which is now exposed as the central star (CS) of the illu- 
minated nebula. Thus the abundances of the nebula should 
reveal information about the chemical processing that took 
place during the AGB; and more precisely, information about 
the last TPs. The spectra of PNe are usually a composite, a 
mixture from the nebula and the illuminating central star. 

Accurate abundances from PNe for helium, C, N, O, Ne, S, 
CI and Ar are available (Kaler 1978; Henry 1989; Dopita et al. 
1997; Stanghellini et al. 2000). These abundances can be used 
as a powerful tool to constrain models of AGB stars (Karakas 
& Lattanzio 2003), in particular the amount of mixing that oc- 
curs during third dredge-up (TDU) episodes following a TP. 
The discovery in PNe of lines of the light neutron-capture el- 
Electronic address: karakas@physics.mcmaster.ca 
Electronic address: M.Lugaro@phys.uu.nl 
Electronic address: gallino@ph.unito.it 

1 Present address: Research School of Astronomy & Astrophysics, Mount 
Stromlo Observatory, Cotter Road Weston Creek, ACT 261 1, Australia 

2 Centre for Stellar & Planetary Astrophysics, Monash University, Clayton 
VIC 3800, Australia 



ement Ge by Sterling, Dinerstein, & Bowers (2002) provides 
not only constraints on the amount of dredge-up but also on 
the operation of the s process occurring in AGB stars. Despite 
the solar Ge abundance being dominated by massive-star nu- 
cleosynthesis, with AGB stars only contributing 6% of solar 
Ge (Arlandini et al. 1999; Busso et al. 2001), this element 
can be produced by the s-process in sufficient quantities to be 
observed in PNe. 

In this Letter we seek to address the question of whether 
detailed models of AGB stars can reproduce the observed 
amount of Ge in PNe. This study is timely given the increas- 
ing number of detections of heavy elements in PNe and in- 
deed this is the first study dedicated to comparing to these 
abundances. In a follow-up paper we will expand upon this 
work to cover all elements up to the first s-peak. 

2. THE NUMERICAL METHOD 

The numerical method we use has been previously de- 
scribed in detail (Lugaro et al. 2004; Karakas et al. 2006). 
Here we summarize the main points relevant for this study. 
We compute the structure first and then perform detailed nu- 
cleosynthesis calculations. The reaction network in the post- 
processing algorithm has been expanded from 74 species to 
156, to include all stable isotopes up to 75 As. We include all 
/3-decays up to arsenic, and all proton, a and n capture reac- 
tions on all species in the network, which results in a total of 
1260 reaction rates. 

To account for neutron captures on elements heavier than 
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TABLE 1 
Details of the stellar models. 



Mass 


Z 


TPs 


7- max 
"ft 


Mass dred 


HBB? 


C/O 


12 C/ 13 C 


Mcnv 


3.0 


0.02 


26 


302 


8.11(-2) 


No 


1.40 


118.0 


0.676 


5.0 


0.02 


24 


352 


5.03(-2) 


Yes 


0.77 


7.840 


1.500 


6.5 


0.02 


40 


368 


4.70(-2) 


Yes 


0.40 


11.60 


1.507 


3.0 


0.012 


22 


307 


9.09(-2) 


No 


2.79 


185.0 


0.805 


6.5 


0.012 


51 


369 


6.48(-2) 


Yes 


0.76 


10.40 


1.389 


1.9 


0.008 


19 


278 


1.70(-2) 


No 


2.08 


143.0 


0.221 


2.5 


0.004 


30 


308 


8.17(-1) 


No 


14.7 


1416 


0.685 



arsenic we use a "double neutron-sink" description (Herwig, 
Langer, & Lugaro 2003), and include two artificial species 
which are linked by the following reactions 75 As(n, 7) 76 g and 
76 g(n,L) 76 g, where 76 g replaces 76 As, and has an initial abun- 
dance equal to the sum of solar abundances from Se to Bi. The 
second artificial particle, L, is equivalent to counting the num- 
ber of neutrons captured beyond arsenic. The ratio (L/ 76 g) is a 
description of the neutrons captured per seed nuclei and could 
in principle be related to the ^-process distribution. 

Most of the 1260 reaction rates are from the 1991 updated 
REACLIB Data Tables (Thielemann, Truran, & Arnould 
1986). Many of the proton, a and n capture-reaction rates 
have been updated by Lugaro et al. (2004). We use the 
NACRE (Angulo et al. 1999) rates for the NeNa and MgAl 
chains, and rates from Karakas et al. (2006) for the 22 Ne +a 
reactions. The cross-section of the 76 g(n,L) 76 g reaction is a 
composite calculated using heavy element distributions pro- 
duced by the s process in low-mass AGB stars from the mod- 
els of Gallino et al. (1998). Test simulations show that using a 
different cross-section for 76 g does not significantly effect the 
elemental abundance of Ge (or any other species) in a 3M Q , 
Z = 0.012 or a 2.5M©, Z = 0.004 model. 

3. THE STELLAR MODELS 

We compute stellar models covering a range of initial mass 
and metallicity listed in the first two columns of Table 1 . The 
metallicities were chosen to reflect the composition of the pro- 
genitor AGB stars which were estimated to range from 0.3Z Q 
to Z Q (Sterling & Dinerstein 2005). The Z = 0.012 stellar 
models were computed with the revised solar elemental abun- 
dances from Asplund, Grevesse, & Sauval (2005) for compar- 
ison to the Z = 0.02 models with abundances from Anders & 
Grevesse (1989). Scaled-solar abundances were assumed for 
the Z = 0.008 and 0.004 models. We use Vassiliadis & Wood 
(1993) mass loss on the AGB. 

In Table 1 we include the number of TPs, the maximum 
temperature in the He-shell, the total mass dredged into the 
envelope during the AGB, the final surface C/O and 12 C/ 13 C 
ratios, and envelope mass at the last computed time-step. All 
masses are in solar units and temperatures in 10 6 K. From Ta- 
ble 1 we see that the 3 and 6.5M Q , Z = 0.012 models were 
similar to their Z = 0.02 counterparts. The low-mass AGB 
models became carbon stars, whereas those models with hot 
bottom burning (HBB) retained an O-rich atmospheric com- 
position (5, 6.5M Q ). There is evidence (e.g. Abia & Isern 
1997) that some extra-mixing takes place in ~ 1 to 3M Q AGB 
stars, resulting in lower observed 12 C/ 13 C and C/O ratios than 
predicted by models such as ours. 

The results presented in this study are subject to many 
model uncertainties including convection and mass loss, 
which affect both the structure and nucleosynthesis. We re- 



fer to Goriely & Mowlavi (2000), Herwig (2005), Ventura & 
D'Antona (2005), and Straniero et al. (2006) for detailed dis- 
cussions on this topic. 

4. THE INCLUSION OF A PARTIAL-MIXING ZONE 

The inclusion of a partial mixing zone (PMZ) at the deepest 
extent of dredge-up will mix protons from the envelope into 
the He-intershell, producing a 13 C pocket. In the PMZ neu- 
trons are liberated during the interpulse period by the reaction 
13 C(a,n) 16 0, and are captured by Fe-seed nuclei to produce 
heavy elements. Observational and theoretical evidence sug- 
gests this is the dominant neutron source in low-mass AGB 
stars (Smith et al. 1987; Gallino et al. 1998). The details 
of how the pocket forms and its extent in mass in the He- 
intershell are still unknown although gravity waves (Denis- 
senkov & Tout 2003), convective overshoot (Herwig 2000), 
induced overshoot during the TDU episodes (Cristallo et al. 
2004; Straniero et al. 2006), and rotationally-induced mixing 
(Langer et al. 1999; Herwig et al. 2003) have been suggested. 

The 22 Ne(a,«) 25 Mg reaction produces a brief strong burst 
of neutrons in the convective pocket during a TP when the 
temperature exceeds ~ 300 x 10 6 K. This is probably the dom- 
inant neutron source in massive AGB stars, where the mass of 
the He-shell is smaller by about an order of magnitude com- 
pared to lower-mass stars, reducing the importance of the 13 C 
pocket. 

As done in previous nucleosynthesis studies (Gallino et al. 
1998; Goriely & Mowlavi 2000; Lugaro et al. 2004) we artifi- 
cially include a PMZ of constant mass. In Table 2 we show the 
maximum mass of the pulse-driven convective region, M CS h, at 
the last TP along with the PMZ masses used. Note that M cs h is 
approximately equal to the He-intershell mass during the TP. 
The PMZ mass is estimated to be between ~ 5% to 15% of the 
mass of the He-rich intershell, and we used PMZs of 1 x 10~ 3 
and 3 x 1O~ 3 M for the 3M Q models. We chose a proton pro- 
file in which the number of protons decreases exponentially 
with the mass depth below the base of the convective enve- 
lope in the same way as described in Lugaro et al. (2004). We 
also computed models with no PMZ for all masses, in order 
to single out the effect of the 22 Ne source. 

5. RESULTS 

In Table 2 we show the surface [Ge/Fe] and [Ga/Fe] ra- 
tios after the last computed TP for each stellar model, where 
we use Fe as the reference element which is justified because 
the surface abundance is predicted by the models to be unal- 
tered by AGB nucleosynthesis 3 . We also provide the surface 
[Ge/S] and [Ge/S] ratios for the 3M , Z = 0.012 model (with 
a PMZ = 2 x 1O~ 3 M0) for comparison. There is no differ- 
ence between these values and those computed using Fe as the 
reference element. Sterling & Dinerstein (2003) tentatively 
detected Ga in one PN (SwSt 1), with the derived [Ga/S] = 
1.64 - 2.87, depending on the level of depletion into dust. Ta- 
ble 2 shows that the largest predicted enhancements in terms 
of [Ge/Fe] were < 0.56 (or a factor of ~ 3.6) for all cases but 
the 2.5M Q , Z = 0.004 model, where [Ge/Fe] = 0.77. Models 
with no 13 C pocket produced very little Ge and Ga, except 
in the case of the 6.5M Q solar metallicity models. For these 
massive AGB models, neutrons are released by efficient acti- 
vation of the 22 Ne neutron source during TPs and the increase 
of Ge at the tip of the TP- AGB is a factor of ~ 1 .5 above so- 
lar. In comparison the 5M , Z = 0.02 model produced little 
Ge with ([Ge/Fe] = 0.1 12) or without (0.04) a PMZ. 

3 Sulfur is similarly not altered. 
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TABLE 2 

[Ge/Fe] and [Ga/Fe] ratios after the last computed TP a 

MASSES ARE IN SOLAR MASSES. 



All 



0.678 0.680 0.682 0.684 0.686 0.688 0.690 0.692 
mass 




0.666 0.668 0.670 0.672 0.674 0.676 0.678 
mass 



FIG. 1. — Composition profile showing the intershell abundances (in log Y, 
where mass fraction X = YA) just before the last computed TP. The shaded 
region is the inner edge of the convective envelope, and a PMZ of 0.002Mq 
was used in both models. We show abundances from (a) the 3M0, Z = 0.012 
model where the Ge intershell abundance after the last TP is 40 Geo, and 
(b) the 2.5Mq, Z = 0.004 model where the Ge intershell abundance is 63 
Ge©. Note that the intershell abundances are typically diluted by one order 
of magnitude at the surface by the last TDU episode. 



The 2.5M Q , Z = 0.004 model produced the most Ge, partly 
because efficient TDU results in a large amount of matter 
dredged from the core into the envelope (see Table 1). An- 
other reason is that the efficiency of the 13 C(a,«) 16 neutron 
source does not depend on the initial Z, and that means more 
neutrons per Fe-seed nuclei are produced at lower metallic - 
ity. Table 1 from Goriely & Mowlavi (2000) also shows more 
Ge in the dredged-up matter with decreasing Z. Gallino et al. 
(1998) Figures 18 and 20 show the s-process enhancement 
factors in the intershell material cumulatively dredged into 
the envelope in 2M Q models with Z = 0.01 and Z = 0.0005, 
respectively. The lower Z model clearly shows an increase in 
the production of the Ge isotopes, with atomic mass ks 70. 

Our results compare favorably with the 3 M Q Z = 
0.02 model computed with the FRANEC evolutionary code 
(Straniero et al. 1997) by Gallino et al. (1998) with their stan- 
dard 13 C pocket (hereafter: the Torino models), who finds 
the final [Ge/Fe], [Ga/Fe], [Zr/Fe] and [Ba/Fe] equal to 0.46, 



Miiss 


x 




PMZ 


[Ge/Fe] 


TGn/Fel 


3.0 


0.02 


1.6(-2) 


2.0(-3) 


0.474 


0.453 


5.0 


0.02 


4.9(-3) 


1.0(-4) 


0.112 


0.111 


6.5 


0.02 


2.3(-3) 


none 


0.191 


0.212 


3.0 


0.012 


1.6(-2) 


1.0(-3) 


0.386 


0.482 








2.0(-3) 


0.557 


0.678 










0.557 b 


0.677 b 


6.5 


0.012 


2.5(-3) 


none 


0.186 


0.304 


1.9 


0.008 


2.0(-2) 


2.0(-3) 


0.252 


0.333 








3.0(-3) 


0.339 


0.439 


2.5 


0.004 


1.7(-2) 


2.0(-3) 


0.777 


0.894 



"Models with no PMZ have [Ge/Fe] es for all cases but the 6.5Mq models. 
b These values are the [Ge/S] and [Ga/S] ratios, computed using sulfur as 
the reference element. 

0.43, 1.03 and 0.93, respectively. The 1.5M Q , Z = 0.02 
produced values equal to 0.28, 0.29, 0.75 and 0.65, respec- 
tively. The Torino 5M Q Z = 0.02 model produced a final 
[Ge/Fe]=0.57, higher than our equivalent model, owing to the 
choice of the Reimers (1975) mass loss (with r\ = 10) which 
results in more TPs (48) compared to our computation (24). 
Massive AGB stars have been suggested as the progenitors 
of Type I bipolar PNe owing to their high He and N/O abun- 
dances (Stanghellini et al. 2006), as well as kinematics (Cor- 
radi & Schwarz 1995). Sterling & Dinerstein (2005) and Ster- 
ling (2006) find a correlation between PNe morphology and 
s-process overabundances, where elliptical PNe are more en- 
riched than bipolar PNe, which show little or no enhancement. 
These observations suggest that Reimers type mass loss is in- 
adequate for modeling massive AGB stars, and formulae with 
a superwind (such as Vassiliadis & Wood 1993) are favored. 

None of the AGB models lost all of their envelopes before 
the computation finished. That is, all of the models were near 
the tip of the AGB but had not proceeded to the post-AGB 
phase. Some of these models could in principle experience 
further TPs and TDUs, and we refer to these as extra TPs. 
The amount of Ge expected after the / th TP can be estimated 
according to X> = (M^X 1 ' 1 + AM Ge )/M* nv , where AM Ge = 
AAMintersheii is the mass of Ge dredged into the envelope, A 
is the TDU efficiency and AMi ntels h e ii is the mass of Ge in 
the intershell. In Figure 1 we show the intershell abundances 
for selected Ge isotopes for two models. A Ge pocket can be 
seen in both cases, and although this will be diluted by one 
order of magnitude by intershell convection, the post-pulse 
abundances are still enhanced compared to solar. 

Using the 3M , Z = 0.012 model with a PMZ of 2 x 
10~ 3 M Q as an example, we estimate that the next TP would 
take place whenM env « 0.7M Q . The mass of the He-intershell 
at the last TP is 0.016M Q , and assuming the same TDU effi- 
ciency as in the previous TP, A = 0.835. At the last computed 
time-step, the Ge envelope and intershell compositions are 
9.64 x 10~ 9 and 1.068 x 10~ 7 (in mass fractions) respectively, 
from which we estimate that after the next TP [Ge/Fe] = 0.64. 
If we assume a mass-loss rate of lO~ 5 M0yr _1 (typical near the 
end of the AGB) multiplied by the interpulse period of 54,000 
years then ~ 0.54M Q is lost per interpulse . From this we 



estimate that the final TP would occur when M„ 



■0.15M 



©• 



We can speculate that if a TDU episode occurs with the same 

4 The same mass model computed with an interpulse period of 25,000 
years near the end of the AGB would allow for 2-3 more TPs. 
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efficiency as before, and using the envelope and intershell Ge 
abundances from the previous estimate (1.17 x 10~ 8 ), we get 
[Ge/Fe] = 0.90 (or Ge = 2.14 x 10" 8 in mass fractions), at the 
tip of the AGB, close to the value required to match the com- 
position of the most Ge-enriched PNe. However caution is 
required because this value is very likely to be an overesti- 
mate, owing to the evidence for a decreasing TDU efficiency 
with a decreasing envelope mass (Straniero et al. 1997). Vari- 
ations in AGB mass-loss rates and model-dependent parame- 
ters such as the interpulse period will also effect this estimate. 

Increasing the size of the PMZ in the models produces more 
Ge but if we increase the size of the pocket too much we run 
into the problem that the PMZ accounts for a large percentage 
of the intershell region. This is contrary to simulations which 
suggest that the 13 C pocket should be at most ~ 10% of the 
intershell (Denissenkov & Tout 2003; Lugaro et al. 2003; Her- 
wig 2005). 

6. DISCUSSION AND CONCLUSIONS 

The results presented in this Letter show that AGB models 
can explain the spread of Ge abundance in PNe from solar to 
a factor of ~ 10 above solar, if we allow for the contribution 
of extra TPs not modeled in detail. Most detailed models pro- 
duced Ge overabundances up to a factor of ~ 3, correspond- 
ing to a dilution of 1 part of He-intershell material per roughly 
20 parts of envelope material. Hence, we can easily account 
for the Ge composition of the majority of the observations, 
in particular taking into account that substantial observational 
uncertainties arise from the choice of the level of dust deple- 
tion. Even for the mostly Ge-enriched PNe (SwSt 1 and BD 
+303639), for which Ge enhancements have been estimated 
to be up to a factor of 30, (but could also be as low as ~ 5 if a 
lower dust depletion is assumed Sterling & Dinerstein 2003, 
2005), it is conceivable that these overabundances could be 
reproduced by considering in more detail the effects of stellar 
model uncertainties, such as the mass loss and treatment of 
mixing processes which affect the TDU in particular. 

The relation between PNe with H-deficient CSs and those 
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with the largest Ge abundances could be explained if a late 
TP during the post-AGB phase (when the envelope mass de- 
creases below ~ 0.01M Q ), was responsible for both. These 
events, described by Werner & Herwig (2006), cause the CS 
to become (partially) H-deficient by mixing and/or burning 
the small remaining envelope with the larger He-intershell, 
which is enriched in He-burning and s-process products. 
However such a scenario would also produce larger enhance- 
ments of the other .s-process elements in this type of PNe, 
while observations show that, on average, the enrichment of 
Kr and Se in PNe with H-deficient CS are not significantly 
different from other PNe (Sterling & Dinerstein 2006; Ster- 
ling 2006). Zhang et al. (2006) recently presented abundances 
for Br, Kr and Xe derived from PNe observations. For all 
of these other elements, which we will discuss in detail in 
a forthcoming paper, results from the Torino code generally 
produce a good match with the observations (except in the 
case of Br which is not observed to be enhanced in PNe con- 
trary to model predictions), with higher enrichments obtained 
at lower metallicities, as in the case of Ge discussed here. 
On the other hand, tentative evidence for the late TP scenario 
could come from BD +303639, with enhanced Ge and an Fe- 
deficient CS (Werner & Herwig 2006), possibly indicating a 
high level of mixing of s-processed material. Further study 
of the stellar model uncertainties and the extra and/or late TP 
scenario together with more accurate observations are needed 
before anything more conclusive can be said. 
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